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In the Minimal Supersymmetric Standard Model (MSSM), the masses of the charged Higgs boson 
(if*) and the CP-odd scalar (A) are related by Mj^ + = M\ + Furthermore, because the 

coupling of W~-A-H + is fixed by gauge interaction, the tree level production rate of qq' —> W ±f —> 
AH ± depends only on one supersymmetry parameter - the mass ( M A ) of A. We show that to a 
good approximation this conclusion also holds at the one-loop level. Consequently, this process can 
be used to distinguish MSSM from its alternatives (such as a general two-Higgs-doublet model) by 
verifying the above mass relation, and to test the prediction of the MSSM on the product of the 
decay branching ratios of A and H ± in terms of only one single parameter - M«. 


PACS numbers: 12.60.-i, 12.15.-y, 11.15.Ex 

One of the commonly discussed new physics models is 
the Minimal Supersymmetric Standard Model (MSSM). 
To describe an experimental data in the framework of the 
MSSM usually requires introducing more than one super- 
symmetry (SUSY) parameters. Hence, the usual prac¬ 
tice is to compare data to a pre-selected class of MSSM 
in which certain well-defined relations among the SUSY 
parameters are assumed in order to reduce the number 
of independent variables needed for discussion. 

An interesting question to ask is “Can one find a pro¬ 
cess to test SUSY models at colliders without making 
many assumptions on the choice of SUSY parameters?” 
To answer that, let us consider the Higgs sector of the 
model. In the MSSM, because of the supersymmetry, 
two Higgs doublets have to be introduced in its Higgs 
sector. Although the MSSM Higgs sector resembles the 
one in a type-II two-Higgs-doublet model (THDM) 0: it 
has a very specific feature - all the Higgs self couplings 
A i are fixed by the electroweak gauge couplings g and g ', 
as required by supersymmetry. Hence, at the tree level, 
only two additional free parameters appear in the Higgs 
sector of the MSSM. We may take Ma (the mass of the 
CP-odd Higgs boson A) and tan (3 (the ratio of the two 
vacuum expectation values) as these two free parameters. 

One of the striking features resulted from the require¬ 
ment of supersymmetry is that the mass (m^) of the 
lightest CP-even Higgs boson (h) has to be less than the 
mass (mz) of the weak gauge boson Z at the Born level, 
although a large radiative correction due to a heavy top 
quark can push this bound up to about 130 GeV in the 
MSSM 0. This result is interesting when compared to 
the theoretical bounds on the mass of the SM Higgs bo¬ 
son. Requiring the SM be a well defined theory up to 
the Planck scale (about 10 19 GeV), the Higgs boson mass 
has to be approximately between 130 GeV and 180 GeV 
00. Therefore, a light Higgs boson with its mass less 
than about 130 GeV can be a signal of the supersym¬ 
metric models, especially the MSSM. However, it is also 
known that such a light Higgs boson can exist in various 
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non-SUSY models, such as a general THDM or the Zee 
model, even when the cutoff scale of the model is close 
to the Planck scale 0. Hence, the existence of a light 
Higgs boson by itself cannot rule out models other than 
the MSSM (or its extensions). 

Another striking feature resulted from the requirement 
of supersymmetry is that the masses of the charged Higgs 
boson H ± and the CP-odd scalar A are strongly corre¬ 
lated. At the Born level, they are related by the mass of 
the W* boson (mw) as 

+ m^y- ( 1 ) 

For comparison, the corresponding mass relation in a 
general THDM is Mj 1± = M\ + |(As — A^v 2 , where v 
is the weak scale (246 GeV) and A^s are two free param¬ 
eters of the model Q. Therefore, the mass relation (jlj) 
can be a strong criterion to discriminate the MSSM from 
its alternatives, e.g., a general THDM. 

To test the mass relation ( 0 ), we propose to study 
the associated production of A and If* at high energy 
hadron colliders, e.g., pp —> AH * at the Fermilab Teva- 
tron (a 2 TeV proton-antiproton collider) and pp —> AH ± 
at the CERN LHC (a 14TeV proton-proton collider). 
As to be explained below, this process has the following 
unique features: (z) its Born level rate generally depends 
on the masses of A and if*. Because of the mass re¬ 
lation ( 0 ), the MSSM prediction of the Born level rate 
only depends on one (in contrast to two or more) SUSY 
parameter - Ma ; (zz) the kinematic acceptance (there¬ 
fore, the detection efficiency) of the signal events do not 
depend on the choice of other SUSY parameters because 
both A and If* are spin-0 (pseudo-)scalar particles so 
that the kinematic distributions of their decay particles 
can be accurately modeled; (in) it can constrain MSSM 
parameters by examining the product of the Higgs boson 
decay branching ratios (in contrast to the product of de¬ 
cay branching ratios and production rate); ( iv ) both M A 
and M h± can be reconstructed from its final state to test 
the mass relation ( 0 ); (v) finally, the electroweak radia¬ 
tive corrections to its production rate and to the mass 
relation ( 0 ) are generally smaller than the expected ex- 


1 



perimental errors, such as the di-jet invariant mass reso¬ 
lution, 

Either in the MSSM or the THDM, the coupling of 
W^-A-H^ is induced from the gauge invariant kinetic 
term of the Higgs sector flj]: 

C in t = %W+(A& i H--H-&*A)+ h.c., (2) 

so that the coupling strength of W + H~ A is completely 
determined by the weak gauge coupling g. (In contrast, 
the coupling constants relevant to the interactions of W- 
boson (or A-boson) and neutral Higgs bosons depend on 
(3 and a, where a characterizes the mixing between the 
two CP-even Higgs bosons h and H.) Thus, the Born 
level production rate of a AH ± pair at hadron colliders 
only depends on M A and M H ±. Since in the MSSM 
these masses are strongly correlated, cf. Eq. ( 0 ), the 
production rate of pp,pp —» AH ± only depends on one 
SUSY parameter, which can be taken as Ma- 

At the Tevatron and the LHC, the dominant con¬ 
stituent process for the production of a AH ± pair is 
qq' —» IU ± * —> AH ± . For a given Ma, the cross section 
cr(pp, pp —*• AH ± ) is completely determined. Its squared 
amplitude, after averaging over the spins and colors is 


\M\ 2 — 7: rn W ( ^FT^ -2 ^2 I-2 p2 

3 (s m - w ) + m w T w 


P 1 sin 2 9 . 


( 3 ) 


where P = \JE\ — M\ with Ea = (s + M\ — 
M^ + )/(2y / s) and 8 is the polar angle of A in the center- 
of-mass (c.m.) frame of A and H ± . We note that for 
the cb —> AH + subprocess, in addition to the CKM 
(Cabbibo-Kobayashi-Maskawa) suppressed s-channel W- 
boson diagram, there is a t-channel diagram, that de¬ 
pends on tan j3. However, the cb —> AH + contribution to 
the inclusive AH + rate is small. For example, its contri¬ 
bution to the total rate is less than 0.01% and 0.1% at the 
Tevatron and the LHC, respectively, for tan (3 = 40 and 
Ma = 90GeV. For a smaller tan (3, its contribution be¬ 
comes negligible. Hence, we shall ignore its contribution 
in the following discussion. In Fig. 1, we show the inclu¬ 
sive production rate of AH ± as a function of Ma- Here, 
the CTEQ5M1 parton distribution functions (PDF) 0] 
are used and both the renormalization and the factor¬ 
ization scales are chosen to be the invariant mass (-y/s) 
of the AH ± pair. A next-to-leading order (NLO) QCD 
correction is also included]^] which typically increases the 
LO rate by about 20% (when the same set of PDF is 
used). 

It is trivial to model the kinematic acceptance (there¬ 
fore, the detection efficiency) of the signal event. This 


1 This is similar to the NLO QCD correction to the W-boson 

production at hadron colliders, except at a different invariant 
mass. 


is because both A and H ± are spin-0 bosons. There¬ 
fore, if the signal is not found, knowing the luminosity of 
the collider, the detection efficiency, and the theoretical 
production rate, one can conclude from the data a con¬ 
straint on the product of the decay branching ratio of A 
and H± as a function of Ma- For example,if the decay 
mode of A — > bb and H + — > t + v t is studied and no ex¬ 
cess is found for a given mass bin of Ma (hence, M H ±) 
when comparing with the experimental data, then one 
can constrain the MSSM by demanding the product of 
the branching ratios, Br(A —> bb) x Br(IL + —> r + ^y), to 
be bounded from above as a function of Ma- Needless to 
say that applying the same strategy, one can constrain 
the product Br(A —> A') x Br(IL + —> Y) for any decay 
mode X and Y predicted by the MSSM as a function of 
only one SUSY parameter Ma- 

In case that a signal is found, the analysis is slightly 
more complicated. In the MSSM, the mass of the heav¬ 
ier CP-even Higgs boson ( H ) is not very different (less 
than about lOGeV) from Ma for Ma ~ 120 GeV and 
tan/3 ~ 10. In this case, qq' —> HH ± can produce the 
similar final states as qq' —» AH ± . Generally, the cou¬ 
pling of W±HH^ depends on g and sin(a — (3). How¬ 
ever, for Ma ~ 190 GeV and tan /3 ~ 10, sin 2 (a — (3) ~ 1 
and the production rate of HH ± is almost the same as 
AH ± in the MSSM. When both of them decay into the 
same decay channels, it will be very difficult to separate 
the production of AH± from HH ± unless a fine mass 
resolution can be achieved experimentally. Nevertheless, 
studying different decay channels can help to separate 
these two production modes. For instance, a heavy H 
can decay into a ZZ pair at the Born level, but A can¬ 
not. 

In conclusion, if no signal is found experimentally, a 
conservative bound on the product of the decay branch¬ 
ing ratios of A and can be derived for a CP- 
conserving model. This is because in a CP-conserving 
model, the AH + and HH + production modes do not in¬ 
terfere even if the masses of A and H are about the same. 
(We note that A is a CP-odd scalar, while H is CP-even.) 

To test the MSSM relation (0) via the AH ± production 
process is straightforward. For example, let us consider 
the bbrv mode. If the signal is sufficiently large as com¬ 
pared to the backgrounds, a resonance bump should be 
observed (at a value (M b i)) in the distribution of the bb 
invariant mass. Then, by searching for the correspond¬ 
ing Jacobian peak (at a value ( Mt{tv ))) in the distribu¬ 
tion of the transverse mass of the tv pair, one can test 
the MSSM by examining whether (. Mt(tv )) is consistent 
with -y/ {M bb ) 2 + within the accuracy of the mass 
resolution of the detector. By testing this mass relation 
via the process pp, pp —> AH± , one can discriminate the 
MSSM from its alternative (e.g., a general THDM). 

In order to prove that the proposed process can be 
used to test the MSSM and is sensitive to only one SUSY 
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FIG. 1. The LO (dotted lines) and NLO QCD (solid lines) 
cross sections of the AH + and AH~ pairs as a function of 
M a . The cross sections for AH + and AH~ coincide at the 
Tevatron for being a pp collider. 


parameter - Ma , we have to show that the SUSY elec- 
troweak correction, which occurs at the one-loop level, 
to the production rate is small. (The final state SUSY 
QCD correction does not contribute until the two-loop 
level.) Specifically, we need to consider two essential 
points: (a) what is the typical size of the radiative cor¬ 
rection to the coupling of W^-A-H^l ( b ) does the mass 
relation (pJ) hold beyond the Born level? In the follow¬ 
ing, we shall have a more detailed discussion on these 
important questions. A brief summary is that as long 
as the typical SUSY mass scale is at the order of a few 
TeV or below, the one-loop correction to the W^AH^ 
vertex can modify the production rate of the AH± pair 
at most by a few percent. Therefore, the electroweak 
correction is smaller than (i) the expected statistical and 
systematic errors of the experimental data, or (ii) the un¬ 
certainty in the theory prediction of the production rate 
originated from the parton distribution functions (which 
is estimated to be about 6% at the Tevatron and 5% at 
the LHC for M A = 120 GeV, when applying the pre¬ 
scription presented in Ref. §), or (iii) the higher order 
(dj or above) QCD correction (which is estimated to be 
about 10% at the Tevatron and less than a percent at 
the LHC for M A = 120 GeV, when varying the factoriza¬ 
tion scale around Vs, the c.m. energy of the subprocess 
qq' —> AiL ± , by a factor of 2). Furthermore, as long as 
the typical SUSY mass scales are at the order of a cou¬ 
ple of TeV or below, the radiative correction to the mass 
relation ( 0 ) is generally smaller than the typical mass res¬ 
olution of the experimental measurement (which is about 
10 GeV for a 100 GeV Higgs boson decaying into jets). 

The dominant one-loop electroweak corrections to the 
qq' —* AH ± process come from the loops of top (t) and 
bottom ( b ) quarks as well as their supersymmetric part¬ 
ners, i.e. stops (< 1 , 2 ) and sbottoms (^ 1 , 2 ), in the MSSM. 
This is due to their potentially large couplings to Higgs 
bosons. 


In the following, we discuss the quark- and squark-loop 
radiative corrections to the effective coupling of W^AH^ 
and to the mass relation (Jl|) . In our calculation, we adopt 
the on-shell renormalization scheme developed by Dabel- 
stein in Ref. [fj (see Appendix A). 

The part of one-loop effective coupling of W^AH^, 
that is relevant to the production process qq' —> H + A , 
can be written asr 


M wha(<?) = ~\(Pa-PhY [l +F (1) ( 9 2 ) 


( 4 ) 


where q V p A and p R are the incoming momenta of W + , 
A and H~, respectively, and g is the effective weak gauge 
coupling evaluated at q 2 . Hence, the radiative correction 
to the cross section of the sub-process qq' —> AH + at the 
one-loop order is 


K^\q 2 ) = 2Re F^{q 2 ). (5) 


The detailed calculation for F^^q 2 ) is summarized in 
Appendix B. As shown in Eqs. (|2l|), ( |49| ) and (|50|), the 
quark-loop contribution is proportional to the squared 
Yukawa coupling constants y 2 (= 2m 2 cot 2 f3/v 2 ) and 
yl(= 2m 2 tan 2 (3/v 2 ). In the large m t or large m;,tan/3 
limit, it can be written as 


pi 1 ) 

quark 


N c 

16-7T 2 


1 


1/3 m 2 1 2 

4» + 2 V S? 1 


( 6 ) 


where N c (= 3) is the number of colors. Since y 2 and y 2 
are at most 0(1) for tan/3 ~ 1 and rrit/mb, respectively, 

^quirk * s mos t a f ew percent for 1 ~ tan/3 ~ mt/rrib. 

We also calculate the squark-loop contribution. As 
compared to the quark effects, the squark effects are 
rather complex due to the additional free (SUSY) pa¬ 
rameters. The mass eigenstates / 1.2 (/ = t or b) of the 
squarks are obtained from the weak eigenstates fp r by 
diagonalizing the mass matrices defined through [^o| 


C 


mass 


~(ft Jr) 


m l m f x f 
m fXf M 2 r 



( 7 ) 


where, M 2 = M? + m 2 + (m| cos2/3)(T/ l —QfSyy) and 
M r = Mq ^ + m 2 + (to| cos2/3)Qfs^y. In this ex¬ 
pression, M?, M?. (for f = t) and M^, (for / = b) 

are the soft-breaking masses for /l, Ir and 1>r, respec¬ 
tively; sw = sin 9w with 0w being the weak mixing an¬ 
gle; Tf and Qf are the isospin and the electric charge 
of the quark /^. Moreover, X t = A t — /xcot/3 and 
Xb = Ab — //tan/3, where A t ( Ab ) is the trilinear A-term 
for t (6), and p, is the SUSY invariant higgsino mass [ |lof . 
For completeness, we have listed all the relevant squark 
and Higgs bosons couplings in Appendix C, so that the 


2 The other form factor, (p A +p H )' x , does not contribute to 
this process for massless quarks. 
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squark-loop contributions to F^^s), cf. (|ii)j), can be di¬ 
rectly read out from the Eqs. (p3|), (|5^ ) and (^4|). 

To examine the effect of one-loop electroweak correc¬ 
tions, we shall discuss two limiting cases below. Firstly, 
we consider the cases with /i = At = A b = 0, i.e., 
the cases without stop mixing (|X t | = 0) and sbottom 
mixing (|V),| = 0). Under this scenario, the masses of 
squarks are proportional to M 2 , and all the relevant cou¬ 
plings between squarks and Higgs bosons are indepen¬ 
dent of the soft-breaking masses Mq, Mg and M g (see 
Appendix C). Thus, the squark-loop effect is decoupled 
and its contribution is very small for a large value of M, 
where M = Mq ~ Mg ~ Mg. (Throughout this pa¬ 
per we denote M as the typical scale of the soft-breaking 
masses.) For a smaller M, F^ ark becomes larger. How¬ 
ever, M cannot be too small because a small M implies 
light squarks whose masses are already bounded from 
below by the direct search results jy|. Furthermore, 
as to be shown later, the case with a small M is also 
strongly constrained by the p parameter measurement. 
Secondly, we examine the case with a large stop mixing, 
assuming m t \ X t ~ M 2 m|. Such a large stop mixing 
leads to a large mass splitting between t\ and ^2 so that 
mg ~ 0(mz) and mg ~ \/2M, while mg 2 ~ M. The 

leading squark contribution to F^^q 2 ) can be expressed 
as 


p(i) 

squark 


-N c 

167T 2 



( 8 ) 


with Y- t = A t cot (3 + p) and = ^p-(Ab tan (3 + p). 
Since in this case |A t | ~ | M 2 /m t ± |p|cot/3|, we have 
\Yi\ ~ 0(M 2 /v) for |/i| ~ M and 1 ~ tan (3. When \ A},\ ~ 
|| > \p\ and tan/3 ~ m t /mb , we find \Y^\ ~ 0(M 2 /v). 
Thus, with a large stop mixing {m t \X t \ ~ M 2 ), A s ( ^ ark 
is proportional to the soft-breaking mass scale M, and 
does not decouple in the large M limit. However, the 
p-parameter (or the T-parameter) can also strongly con¬ 
strain such kind of model. With a large stop mix¬ 
ing (M 2 ~ m t \X t \), the squark contribution (cf. Ap¬ 
pendix D) to the p-parameter is 


Ap squark ~ (2.2 x 10 -3 )— 2 ". 


(9) 


Since any new physics contribution to the p-parameter 
has to be bounded by data as (jTl|] 

— 1.7 < Ap n ew x 10 3 < 2.7, at 2 ct level, (10) 


the scale M cannot be too large in this case. Conse¬ 
quently, the above F s ^ ark is constrained to be smaller 
than a few percent as long as p 2 is not much larger than 
M 2 . 

To examine the effect from the stop and sbottom loops 
to the production rate of AH ± , we consider 4 sets of 
SUSY parameters, as listed in Table 1, which give the 



Setl 

Set2 

Set 3 

Set4 

Mq =M 0 = M b (GeV) 

106 

84 

408 

409 

tan f3 

2 

40 

2 

40 

At = A b (GeV) 

0 

0 

+1261 

+1119 

H 

0 

0 

+300 

+300 

mg (GeV) 

197 

184 

100 

100 

mg (GeV) 

199 

188 

612 

611 

mg (GeV) 

108 

88 

407 

373 

mg (GeV) 

116 

103 

412 

447 

^Psquarks X 10 

2.72 

2.70 

2.71 

2.70 


TABLE I. The SUSY (input and output) parameters used 
in Fig. 2. 


largest allowed deviation in the p-parameter. Set 1 and 
Set 2 represent the cases without either a stop mix¬ 
ing ( X t = 0) or a sbottom mixing (Xb = 0), and 
Set 3 and Set 4 are the cases with a large stop mixing 
(m t \X t \ ~ M 2 ) and mg ~ 100 GeV. The A'W(s) factor, 
as defined in Ecp ([)]), is shown in Fig. 2 as a function 
of the invariant mass (-y fs ) of the constituent process for 
M a = 90 GeV. It is clear that the quark-loop contribu¬ 
tion to /\H)(g) dominates the squark-loop contribution. 
For the above sets of SUSY parameters (Set 1-4), the 
squark-loop contribution is smaller than the quark-loop 
contribution by about a factor of 100. Generally, the 
squark contributions are at most a few percent, unless 
|p| is taken to be very large as compared to the scale M. 
We have checked that this conclusion does not change 
when our assumption of Mq ~ M^ ~ M^ is relaxed to 
some extent. Including both the quark- and squark-loop 
contributions to A'G^(s), we found that the correction to 
the hadronic cross section of H + A production in the in¬ 
variant mass region just above the H + A threshold, where 
the constituent cross section is the largest, is at a percent 
level. In summary, we illustrated that to be consistent 
with the low-energy data and the direct search results 
for stops and sbottoms, the one-loop electroweak correc¬ 
tion to the production rate of pp, pp —> AH ± is small (at 
most a few percent). 

Next, we discuss the one-loop corrections to the mass 
relation 0) Let us parameterize the deviation from the 
tree-level relation by 5 , so that at the one-loop order 

M h ± = ^Jm\ + m 2 w (i + S). (11) 

We note that in our renormalization scheme (see Ap¬ 
pendix A), Ma and mw are the input parameters, but 
M h + is not. The one-loop corrected mass of the charged 
Higgs boson M H± can be obtained by solving 


0 = Det 


4V(p 2 ) 4V(p 2 ) 

t^ g _(t 2 ) r %\ h Ap 2 ) 


( 12 ) 


where T,- 2 ^(p 2 ) represent the renormalized two-point 
functions in the basis of the renormalized Goldstone bo¬ 
son (G ± ) and charged Higgs boson (H ± ) fields. Here, 
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FIG. 2. The A-factor, A ( 1 ) (s), of the constituent process 
qq' — > H + A for M A = 90 GeV, as a function of the invariant 
mass y/s of qq . The solid lines correspond to the top and 
bottom quark contribution. The cases where the squark-loop 
contribution is included are described by the dotted lines for 
those without stop mixing (Set 1 and Set 2) and by dashed 
lines for those with maximal stop mixing (Set 3 and Set 4), 
respectively. 


The notation “Det” denotes taking the determinant of 
the 2x2 matrix. One of the solution of the above equa¬ 
tion is p 2 = 0, which corresponds to the charged Nambu- 
Goldstone mode, and another is At the one-loop 

level, the pole mass of the charged-Higgs boson can be 
calculated from 

Mfl± = M A + niw 

-\~n AA (M A ') ~ n h+H~ (M\ -f- 1Tl\y) + Tlww( m w)i (13) 

where II<^(g 2 ) {<j> = A,H ± ,and W) are the self-energies. 
For completeness, we list the quark and squark contribu¬ 
tions to the self-energies of A and in Appendix E. 

When Atfi and p are zero (i.e., no-mixing case), the 
leading contribution (which is proportional to the forth 
power of heavy quark mass) to S is found to be 

s ~ N c f m t m l \ 1 

8ir 2 v 2 \ M\ + m^ J sin 2 /? cos 2 f3 

This correction is substantial for tan (3 ~ mt/mb and 
M 2 m 2 . Applying Eq. © with the complete ex¬ 
pression of II ^(q 2 ), we found S to be less than 4.9% for 
2 < tan/3 < 40, M < 2000GeV and M A > 90 GeV. Our 
result agrees well with Ref. ju|, in which the approxi¬ 
mate formula were presented for M 2 m 2 . 

For the cases with a nonzero Ab,t and p, S receives 
extra contributions, which are proportional to A 3 4 b /M 4 , 
A 2 b p 2 /M 4 and p 4 /M 4 Jfj.O] originated from the squark 
couplings [cf. Eqs. (p9|), (pOf), ( |33| ) and ([m])] and squark 
masses. For the Set 3 and Set 4 parameters listed in Ta¬ 
ble 1, S is less than 5.3% and 3.6% for M A > 90 GeV, 
respectively. In summary, as long as and |/i| are 



not too large as compared to M, in a wide range of the 
parameter space that is allowed by the available exper¬ 
imental and theoretical constraints, <5 does not exceed 
7-10%. 

Supported by the finding that the one-loop electroweak 
corrections to the W^AH^ coupling and to the mass 
relation M‘j I± = M\ + rriyy are generally smaller than 
the other theoretical errors (such as the parton distri¬ 
bution function uncertainties) and the expected experi¬ 
mental errors (such as the mass resolution of Higgs bo¬ 
son decaying into jets), we anticipate that our conclu¬ 
sions based upon a Born level analysis should also hold 
well at the loop level. Namely, studying the process 
pp,pp —> W ±x —> A7T ± allows us to distinguish the 
MSSM from its alternatives by verifying the mass re¬ 
lation (fl;) and checking its production rate. If a signal 
is not found, studying this process can provide an upper 
bound on the product of the decay branching ratios of A 
and IT ± as a function of the only one SUSY parameter 
- M A . 

To detect the signal event, it is necessary to suppress 
its potentially large backgrounds. For example, the rvbb 
backgrounds can be largely reduced by having a good 
6-tagging and tau selection (by using the nature of r po¬ 
larization, which differs between a parent H ± and IV ± 
0 ). We expected that its observability is relatively easy 
at the LHC and is a challenging task at the Tevatron for 
the signal event rate is small. Clearly, a detailed Monte 
Carlo analysis is needed to calculate the significance of 
the signal event at a collider. This will be deferred to a 
future study. 


We thank K. Hagiwara, G.-L. Kane, P. Nadolsky, 
Y. Okada and T. Tait for useful discussions. This work 
was supported in part by the NSF grant PHY-9802564. 

APPENDICES 

A. Renormalization 

In this paper, we adopt the on-shell renormalization 
scheme developed by Dabelstein [|J to calculate the one- 
loop electroweak corrections. The standard model pa¬ 
rameters are fixed by defining a em , rtiw and mz, and 
the additional SUSY parameters in the Higgs sectoi^jare 


3 There are 7 parameters in the Higgs sector of the MSSM. 
They are g ', g , vi, t> 2 , mi, m 2 , and m 3 . Beyond the Born 
level, the wavefunction renormalization factors Zhi and Zh 2 
also need to be introduced to renormalize the theory, where 

Hi and H 2 denote the two Higgs doublets in the model. 
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fixed by the following renormalization conditions: (1) the 
tadpole contributions (Th 1 = 0, Th 2 = 0), (2) the on- 
shell condition for the mass of A , (3) the on-shell con¬ 
dition for the wavefunction of A , (4) a renormalization 
condition on tan/3 (which requires 6v\/v\ = 5v 2 /u 2 ), and 
(5) a vanishing A — Z mixing for an on-shell A. 


B. Calculation of F < ' 1 \q 2 ) 


The one-loop correction to the renormalized form fac¬ 
tor of the W ± H^ L A vertex, apart from the effective weak 
gauge coupling g , can be written as 

F^(q 2 ) = ~zTa^I h - {1 + SF WHA 

+Fwha{M a , Mh, g 2 )} — 1, (15) 


where Z AA and Z H + H - are the finite wavefunction fac¬ 
tors for the renormalizations of the external Higgs bosons 
A and H±. In our scheme, 

Z AA = 1, (16) 

Z H+H - = 1 - W H+H _ (Mi + m 2 w ) + W AA (M 2 A ) , (17) 


where H' AA (M A ) denotes taking the derivative of the two 
point function Il AA (k 2 ) of the CP-odd scalar A with re¬ 
spect to k 2 at k 2 = M 2 , etc. The terms inside the 
curly bracket of Eq. ( |l5| ) arise from the renormalized 
vertex function of WFIA. F^jj A (p 2 A ,p 2 H ,q 2 ) represents 
the one-loop contribution of the one-particle-irreducible 
(1PI) diagrams with p 2 A ,p 2 Hl q 2 as the four-momentum 
square of the incoming A, and particles, respec¬ 
tively. SFwha is the counterterm contribution resulting 
from the field renormalization of F[ + and A: 


H+A -> H+A 


1 + -SZ H + + -SZ A 


(18) 


In terms of the independent counterterms fixed by 
the renormalization scheme, the wavefunction countert¬ 
erms SZ H + and 5Z A can be written as (sin 2 (3)SZh x + 
(cos 2 (3)5Z H2 which is found to be equal to — |n AA (M A ). 
We note that in SFwha the contributions from the coun¬ 
terterms of the weak gauge coupling and the wavefunc¬ 
tion renormalization of the W-boson are not included, 
because they should be combined with the W-boson self 
energy contribution to derive the running weak gauge 
coupling g(q 2 )- In our numerical calculation, we use 

g 2 = 4V2mwG F . 


In summary, the one-loop electroweak correction to 
F^V) is found to be 

F^(q 2 ) = F^ A (M 2 A ,M 2 H± ,q 2 ) 


~2^' h+h ~ (^a + «■%) — 2 U> aa( m a)- ( 19 ) 


In the above equation, the top- and bottom-loop contri¬ 
bution to FyyHA I s gi yen by 

= *52 Fwha^’PA’Ph), ( 20 ) 

f f f'=ttb,bbt 

with 

pi if „2 2 \ _ , N C 2 f 2 r fif 2 r fff 

PwhaVPaiPhiQ ) - \P A ^31 ~ Ph^ 32 

+(2 PA ■ PH - P A )Cll f - {2p A ■ p H - p 2 H )Cil f 
+(D + 2- Cii f ') +p 2 A C({ f ' - (2p A -ph 
+Ph)cI l f ' - 2p\cH f ' -(D- 2)C 2 / / / ' - m 2 C{[ s ' 

- ( q 2 + to 2 ) C{l f } - Cfj^y f yf'rn f mf,Cf ff . (21) 

where c/ = +1 and —1 for ///' = ttb and bbt , respec¬ 
tively, and C{j* are defined in terms of the Passarino- 
Veltman functions [|l5| with 

C 'i/ f = Ci i (p%Ph> (pa + Ph) 2 ', m/, to/, to//) . (22) 

The stop- and sbottom-loop contribution is given by 

2 

plPI(squark), 2 2 2 \ _ N c 

+WHA \PAiPHi 0 )— + l 

i,j } k=l 

X {u* L D Lk (i\[tj,ti,A]) X[b* k ,tj,H~}C iiijbk 

- U* kL D Li ( iX[b *, bj,A]) X[b*,t k ,H-}C lil f k } , (23) 

where Un, Du are the rotation matrices for stops and 
sbottoms between the weak eigenstate basis (I = L, R) 
and the mass eigenstate basis (i = 1,2), respectively. 
A[/*,/), 4>k] represents the coefficient of the f*fj<f>k inter¬ 
action in the MSSM Lagrangian, as listed in Appendix C, 
and 

(jiii'jik = (c n - C 12 ) (pi,P2,<?‘,rnj i ,rrijpmj„') . (24) 

The quark (top and bottom) and squark (stop and sbot- 
tom) loop contributions to the self-energies T3aa(<1 2 ) and 
H h+h- (q 2 ) can be found in Appendix E. 


C. Squark couplings with H+ and A 


The mass eigenstates of the squarks are related to their 
weak eigenstates by the rotation matrix Oj) with R = 
Ei 0{^i f i, where i = 1,2 and I = L,R; = Un and 
Dj i for / = t and 6, respectively. In terms of the mixing 
angles Of, we have 


fi \ = ( cos 0 f sin Of \ ( f L \ 

h) V “ sin0 / cos6 f ) V h / 


(25) 


Here, we define the mixing angle Of so that fi is lighter 
than / 2 . 
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The coupling constants among the weak-eigenstate 
squarks and the Higgs bosons are defined through the 
Lagrangian 




(26) 


D. Squark contributions to the p parameter 

The squark one-loop contribution to the p parameter 
is given by 

Ap = p- 1 = -4v^G F Re[An^(0) - An| 3 (0)], (45) 


Hence, the coupling constants for the mass-eigenstate 
squarks are a linear combination of the couplings for the 
weak-eigenstate squarks, and 

a [/?,/', -] = al oQfj , 0 ,...] 

= O f u O{'J\[fI J'j, (27) 

The relevant couplings X[fj , fj, </>,...], denoted as 
Ai ar e listed below. 


with 117 


AnW) = is?£ E T /dioLi 2 |ofii 2 

f=t.,bi,j =1 


xS(g 2 ;m|,mp, (46) 


AH 33 ^) =-E_ |b r Li| 2 |Il{ J | 2 H((7 2 ;m|.,m| ), (47) 


32tt 2 




A 


—\/2 

- (m^ sin 2/3 — to 2 tan /3 — m 2 cot /3), (28) 


*>3=i 

where O),- = Un and O b H = D/*; B{q 2 - 1 m\,ml) = 
H(m 2 ) + A(m|) — 3B 22 {q' 2 -, to 2 ,m|). By using the ex¬ 
pression 

to, l,o 0, m?77io , mo 

B(0;m?,m^) =--(m? + m|)-l- 2 1 2 ln ~I . ( 4S ) 

2 777-f — 777-2 m i 

Eq. © is deduced under the assumption that M 2 = 
Mq ~ My ~ M 2 to 2 and the stop mixing is large 
(?7Ti|X t | ~ M 2 and m.b|Xb| ~ 0), so that ~ 0{mz), 
which yields mj 2 ~ V2 M, and m bi ~ mj ~ M. 


v 

V2m t 


{A t cot/3-f /z), 


b* L tRH y 

. v2 m b , . _ , 

A i* R t L H- = —^—( Ab tan /3 + W, 

. 2 \[2m t m b 

b* R t R H- ~ v g in 2 /3 ’ 

= A/*/ h A = °> (/ = *! A 

= — (H b tan/3 + /i), 

. imt , . „ x 

A** = —I 24 * cot P + W. 


t; 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


2 2 

hib L AA = -^ptan 2 /3+ ^-(I), - QbSw) cos 2/3, (35) 

2 2 

AgjS*^ = tan2 P + Ef C0S 2/? ’ 

2 2 

X bii L AA = —^T cot2 P+^f( T t ~ Qt s w) cos 2/3, (37) 

2 9 

_ 777 ,^ (7^* 

A^b B AA = cot2 P + y cos 2/3, 


(36) 


(38) 


E. Self energies 

The (top and bottom) quark-loop contributions to the 
self-energies n^^g 2 ) and A H + H -{q 2 ) are expressed in 
terms of the Passarino-Veltman functions [^5) by 

^V) = E 2 V 2 f U 2 {B 1 {q 2 ,m f ,m f ) 


f=t,b 


A 


b* T b L H+H- 


_O 770 2 Cl 2 

*■ tan 2 /3 + yA(T t — Qts^) cos 2/3, (39) 


Ar 


2 

2?T7.? 

tan 2 (3 + ErQbS^ cos 2/3, (40) 

_2?77 ^ 

t-j L H+H- = yE tan2 P + 2 ^ cos2/? ’ ( 41 ) 


b* R b R H+H- - y 2 UCU1 e>-r 


+B 21 {q 2 ,m f ,m f )) + L> • B 22 {q 2 , m f , m f ) 
+m 2 f B 0 {q 2 ,m f ,m f )} , (49) 

H^Th- (l 2 ) = - yy 2 + ?/*) {? 2 (-Sife 2 . mb, m t ) 

+B 2 i{q 2 ,m. b ,m t )) + D ■ B 22 {q 2 , m b , m t )} , 

N r 


167T 2 


4y b y t m b m t B 0 {q , m b , m t ). 


(50) 


_ 2777 ,^ 

i* R i R H+H- = cot 2 /3 + -Z-QtSw cos 2/3, 


(42) 


The stop- and sbottom-loop contributions are given by 

iV, 2 


n'trV) = -^5 E E 


where T t ,T b , Qt and Qb are A, -A, | and -A, respectively, 
and 


f=t,b *>3=i 
,2 


X t}bjH+ — As*tjff-> 
= ~ x r L f R Ai 


(43) 

( 44 ) 


xA[/*,/j,H]A[/*,/i,H]H 0 te (51) 

iV, 


16-7T 2 


2 E E^./bA^m^), 


( 52 ) 


/=t,b *=i 


ns;"( ? 2 ) = -A E 


for I,J = L 1 R and f = t,b. 


167T 2 


*>3=i 
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(53) 


xA [t *, bj , H + }X[b* , ti , H ]B 0 {q 2 ,mi.,m- b .) 

N . 2 

-i( 54 ) 

f=t,b * =1 

The self-energy II ww {q 2 ) of the IT boson was already 
presented in the literature. For example, the quark-loop 
contribution can be found in Ref. [(Tb), and the squark- 
loop contribution in Ref. []l7| . 
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